Archives of Toxicology
○ Springer Science and Business Media LLC
All preprints, ranked by how well they match Archives of Toxicology's content profile, based on 18 papers previously published here. The average preprint has a 0.01% match score for this journal, so anything above that is already an above-average fit. Older preprints may already have been published elsewhere.
Wills, J. W.; White, A.; Harte, D. S. G.; Buckley, R.; Harvey, J. S.; Lynch, A. M.
Show abstract
The management of N-nitrosamine impurities challenges pharmaceutical development and regulation worldwide. Because most medicinal exposures are shorter than lifetime and absolute impurity exclusion is impossible, reliable approaches to define duration-specific intake limits are essential. On the premise that carcinogenic risk is proportional to cumulative dose, the Less-Than-Lifetime (LTL) Threshold of Toxicological Concern (TTC) framework defines progressively lower intake limits for mutagenic impurities over longer exposures. However, N-nitrosamines are currently treated as a cohort of concern, necessitating compound-specific evaluation placing reliance on in vivo mutagenicity assays for impurity qualifications. To better understand durational potency relationships and the application domain of the LTL-TTC, we apply benchmark dose (BMD) modelling to cumulative-dose-scaled transgenic rodent (TGR), error-corrected sequencing and rodent carcinogenicity datasets for N-nitrosodimethylamine (NDMA) obtained from the published literature. For TGR, cumulative-dose scaling better resolved liver as the most sensitive organ and reduced interstudy variability: liver BMDs spanned [~]80-fold in daily-dose units but only [~]20-fold when scaled to cumulative dose. Among closely-matched mouse liver gavage studies, cumulative-dose BMDs only varied by [~]2.5-fold across 1 to 28-day treatment regimens. Error-corrected sequencing also demonstrated parity, with acute-dose regimens producing mutation burdens near-identical (< 1.2-fold) to those cumulated from 28-day repeat-dose regimens. Comparable results were obtained from carcinogenicity datasets confirming proportionality-of-effect to cumulative dose. These findings empirically support the validity of the LTL-TTC concept. More broadly, they demonstrate that short-term in vivo mutagenicity assays can serve as reliable surrogates for lifetime carcinogenicity studies, strengthening the scientific and regulatory basis for duration-adjusted acceptable intakes for N-nitrosamine impurities.
Edizcan, Z.; Schaller, S.; Kuepfer, L.; Geci, R.
Show abstract
High-throughput physiologically based kinetic (HT-PBK) modelling provides a mechanistic framework for predicting systemic pharmacokinetics (PK) from in vitro and in silico data, supporting non-animal chemical safety assessment within Next generation risk assessment (NGRA). Here, we applied HT-PBK modelling to dermal exposure, a key route of human contact with chemicals. Using the skin permeation model of the Open Systems Pharmacology Suite, we simulated systemic PK profiles of 52 compounds based solely on physicochemical properties predicted by quantitative structure-activity relationship (QSAR) models and without any compound-specific in vitro measurements. We systematically compared different QSAR tools for lipophilicity, solubility, and other parameters to identify optimal model parameterisation strategies. Across all compounds, the best-performing HT-PBK strategy predicted 75% of Cmax and AUC values within a tenfold range of observed human plasma data extracted from published clinical studies. A systematic tendency toward overprediction of systemic PK was observed, likely due to missing study metadata and the default assumption of fully hydrated skin. Prediction errors were larger for dermal than for oral exposure, reflecting the greater complexity and variability of dermal absorption processes. Nevertheless, key exposure metrics were reproduced within acceptable limits. These results demonstrate the feasibility of fully in silico, non-animal HT-PBK modelling for dermal absorption and support its use as a pragmatic tool for exposure-driven safety assessment within NGRA frameworks.
Niemeijer, M.; Wiecek, W.; Huppelschoten, S.; Bouwman, P.; Baze, A.; Parmentier, C.; Richert, L.; Paules, R. S.; Bois, F. Y.; van de Water, B.
Show abstract
Background & AimsOne of the early key events of drug-induced liver injury (DILI) is the activation of adaptive stress responses, a cellular mechanism to overcome stress. Given the diversity of DILI outcomes and lack in understanding of population variability, we mapped the inter-individual variability in stress response activation to improve DILI prediction. Approach & ResultsHigh-throughput transcriptome analysis of over 8,000 samples was performed in primary human hepatocytes of 50 individuals upon 8 to 24 h exposure to broad concentration ranges of stress inducers: tunicamycin to induce the unfolded protein response (UPR), diethyl maleate for the oxidative stress response, cisplatin for the DNA damage response and TNF for NF-{kappa}B signalling. This allowed investigation of the inter-individual variability in concentration-dependent stress response activation, where the average of benchmark concentrations (BMCs) had a maximum difference of 864, 13, 13 and 259-fold between different hepatocytes for UPR, oxidative stress, DNA damage and NF-{kappa}B signalling-related genes, respectively. Hepatocytes from patients with liver disease resulted in less stress response activation. Using a population mixed-effect framework, the distribution of BMCs and maximum fold change were modelled, allowing simulation of smaller or larger PHH panel sizes. Small panel sizes systematically under-estimated the variance and resulted in low probabilities in estimating the correct variance for the human population. Moreover, estimated toxicodynamic variability factors were up to 2-fold higher than the standard uncertainty factor of 101/2 to account for population variability during risk assessment, exemplifying the need of data-driven variability factors. ConclusionsOverall, by combining high-throughput transcriptome analysis and population modelling, improved understanding of variability in stress response activation across the human population could be established, thereby contributing towards improved prediction of DILI.
Hardy, B.; Mohoric, T.; Parween, S.; Lopez, B.; Daligaux, P.; Darde, T.; Chesne, C.; Stockman, N.; Lemos, J.; Saiakhov, R.; Boglari, C.; Poon, A.; Ukaegbu, D.; Andrade, C.
Show abstract
The transition toward animal-free safety assessment of chemicals has accelerated the development of New Approach Methodologies (NAMs) for predicting skin sensitization. However, individual in silico models and experimental NAM assays frequently produce inconsistent or contradictory results, limiting their reliability when used in isolation. To address this challenge, we present a tiered integrated assessment framework implemented through the open source SaferSkin application, which enables systematic comparison and integration of multiple predictive models and experimental data within a transparent weight-of-evidence workflow. In this case study, a diverse set of 21 reference compounds was evaluated using a battery of in silico approaches, including the OECD QSAR Toolbox, VEGA, CASE Ultra and additional machine-learning models implemented within SaferSkin. The platform enables side-by-side comparison of predictions and integration of experimental data through Bayesian network models, allowing probabilistic updating of predictions as new evidence becomes available. Our results demonstrate that reliance on any single predictive model is insufficient for robust hazard identification due to frequent disagreement between models. In contrast, consensus interpretation across multiple modelling approaches combined with targeted experimental evidence substantially improves predictive confidence. The integrated weight-of-evidence framework showed strong concordance with reference classifications and was further supported by independent validation using the Pred-Skin Bayesian model. Importantly, the tiered workflow enables resolution of ambiguous cases. For example, lower-tier predictions for ethyl (2E,4Z)-deca-dienoate were inconsistent across models, whereas targeted third-tier testing using the SENS-IS assay identified the compound as a strong sensitiser (GHS Category 1A). Overall, this study demonstrates how integrated modelling, Bayesian evidence updating and targeted NAM testing can reduce uncertainty in skin sensitization assessment. The SaferSkin framework provides a transparent and reproducible approach for implementing Next Generation Risk Assessment (NGRA) strategies and supports the development of animal-free regulatory toxicology and Safe-and-Sustainable-by-Design chemical innovation. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=115 SRC="FIGDIR/small/711911v1_ufig1.gif" ALT="Figure 1"> View larger version (45K): org.highwire.dtl.DTLVardef@b59ca0org.highwire.dtl.DTLVardef@13de455org.highwire.dtl.DTLVardef@599358org.highwire.dtl.DTLVardef@d87fd1_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOGraphical AbstractC_FLOATNO C_FIG
Rowan-Carroll, A.; Readon, A.; Leingartner, K.; Gagne, R.; Williams, A.; Kuo, B.; Bourdon-lacombe, J.; Moffat, I.; Carrier, R.; Nong, A.; lorusso, l.; Ferguson, S. S.; Atlas, E.; Yauk, C.
Show abstract
Per- and poly-fluoroalkyl substances (PFAS) are widely found in the environment because of their extensive use and persistence. Although several PFAS are well studied, most lack toxicity data to inform human health hazard and risk assessment. This study focussed on four model PFAS: perfluorooctanoic acid (PFOA; 8 carbon), perfluorobutane sulfonate (PFBS; 4 carbon), perfluorooctane sulfonate (PFOS; 8 carbon), and perfluorodecane sulfonate (PFDS; 10 carbon). Human primary liver cell spheroids (pooled from 10 donors) were exposed to 10 concentrations of each PFAS and analyzed at four time-points. The approach aimed to: (1) identify gene expression changes mediated by the PFAS; (2) identify similarities in biological responses; (3) compare PFAS potency through benchmark concentration analysis; and (4) derive bioactivity exposure ratios (ratio of the concentration at which biological responses occur, relative to daily human exposure). All PFAS induced transcriptional changes in cholesterol biosynthesis and lipid metabolism pathways, and predicted PPAR activation. PFOS exhibited the most transcriptional activity and had a highly similar gene expression profile to PFDS. PFBS induced the least transcriptional changes and the highest benchmark concentration (i.e., was the least potent). The data indicate that these PFAS may have common molecular targets and toxicities, but that PFOS and PFDS are the most similar. The transcriptomic bioactivity exposure ratios derived here for PFOA and PFOS were comparable to those derived using rodent apical endpoints in risk assessments. These data provide a baseline level of toxicity for comparison with other known PFAS using this testing strategy.
Roe, H.; Tsai, H.-H. D.; Ball, N.; Wright, F. A.; Chiu, W.; Rusyn, I.
Show abstract
An important element of the European Unions "Registration, Evaluation, Authorisation and Restriction of Chemicals" (REACH) regulation is the evaluation by the European Chemicals Agency (ECHA) of testing proposals submitted by the registrants to address data gaps in standard REACH information requirements. The registrants may propose adaptations, and ECHA evaluates the reasoning and issues a written decision. Read-across is a common adaptation type, yet it is widely assumed that ECHA often does not agree that the justifications are adequate to waive standard testing requirements. From 2008 to August 2023, a total of 2,630 Testing Proposals were submitted to ECHA; of these, 1,538 had published decisions that were systematically evaluated in this study. Each document was manually reviewed, and information extracted for further analyses. Read-across hypotheses were standardized into 17 assessment elements (AEs); each submission was classified as to the AEs relied upon by the registrants and by ECHA. Data was analyzed for patterns and associations. Testing Proposal Evaluations (TPEs) with adaptations comprised 23% (353) of the total; analogue (168) or group (136) read-across adaptations were most common. Of 304 read-across-containing TPEs, 49% were accepted; the odds of acceptance were significantly greater for group read-across submissions. The data was analyzed by Annex (i.e., tonnage), test guideline study, read-across hypothesis AEs, as well as target and source substance types and their structural similarity. While most ECHA decisions with both positive and negative decisions on whether the proposed read-across was adequate were context-specific, a number of significant associations were identified that influence the odds of acceptance. Overall, this analysis provides an unbiased overview of 15 years of experience with testing proposal-specific read-across adaptations by both registrants and ECHA. These data will inform future submissions as they identify most critical AEs to increase the odds of read-across acceptance.
Reardon, A. J. F.; Rowan-Carroll, A.; Ferguson, S. S.; Leingartner, K.; Gagne, R.; Kuo, B.; Williams, A.; Lorusso, L.; Bourdon-Lacombe, J.; Carrier, R.; Moffat, I.; Yauk, C. L.; Atlas, E.
Show abstract
Per- and polyfluoroalkyl substances (PFAS) are some of the most prominent organic contaminants in human blood. Although the toxicological implications from human exposure to perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA) are well established, data on lesser-understood PFAS are limited. New approach methodologies (NAMs) that apply bioinformatic tools to high-throughput data are being increasingly considered to inform risk assessment for data-poor chemicals. The aim of this investigation was to identify biological response potencies (i.e., benchmark concentrations: BMCs) following PFAS exposures to inform read-across for risk assessment of data-poor PFAS. Gene expression changes were measured in primary human liver cell microtissues (i.e., 3D spheroids) after 1-day and 10-day exposures to increasing concentrations of 23 PFAS. The cells were treated with four subgroups of PFAS: carboxylates (PFCAs), sulfonates (PFSAs), fluorotelomers, and sulfonamides. An established pipeline to identify differentially expressed genes and transcriptomic BMCs was applied. We found that both PFCAs and PFSAs exhibited a trend toward increased transcriptional changes with carbon chain-length. Specifically, longer-chain compounds (7 to 10 carbons) were more likely to induce changes in gene expression, and have lower transcriptional BMCs. The combined high-throughput transcriptomic and bioinformatic analyses supports the capability of NAMs to efficiently assess the effects of PFAS in liver microtissues. The data enable potency ranking of PFAS for human liver cell spheroid cytotoxicity and transcriptional changes, and assessment of in vitro transcriptomic points of departure. These data improve our understanding of the health effects of PFAS and will be used to inform read-across for human health risk assessment.
Cho, E.; Swartz, C. D.; Williams, A.; Rivas, M.; Recio, L.; Witt, K. L.; Schmidt, E. K.; Yaplee, J.; Smith, T.; Van, P.; Lo, F. Y.; Valentine, C. C.; Salk, J. J.; Marchetti, F.; Smith-Roe, S.; Yauk, C.
Show abstract
Error-corrected Duplex Sequencing (DuplexSeq) enables direct quantification of low-frequency mutations and offers tremendous potential for chemical mutagenicity assessment. We investigated the utility of DuplexSeq to quantify induced mutation frequency (MF) and spectrum in human lymphoblastoid TK6 cells exposed to a prototypical DNA alkylating agent, N-ethyl-N-nitrosourea (ENU). Furthermore, we explored appropriate experimental parameters for this application, and assessed inter-laboratory reproducibility. In two independent experiments in two laboratories, TK6 cells were exposed to ENU (25-200 {micro}M) and DNA was sequenced 48, 72, and 96 h post-exposure. A DuplexSeq mutagenicity panel targeting twenty 2.4-kb regions distributed across the genome was used to sample diverse, genome-representative sequence contexts. A robust increase in MF that was unaffected by time was observed in both laboratories. Concentration-response in the MF from the two laboratories was strongly positively correlated (R2=0.95). C:G>T:A, T:A>C:G, T:A>A:T, and T:A>G:C mutations increased in consistent, concentration-dependent manners in both laboratories, with high proportions of C:G>T:A at all time points. The target sites responded similarly between the two laboratories and revealed a higher average MF in intergenic regions. These results, demonstrating remarkable reproducibility across time and laboratory for both MF and spectrum, support the high value of DuplexSeq for characterizing chemical mutagenicity in both research and regulatory evaluation.
Lomash, V.; Srinivasan, M.; Pitthala, M.; Sayeed, A.; Venkatesan, G.; Joseph, B.
Show abstract
Evaluation of unintended immunotoxicity represents an important component of nonclinical safety assessment, as perturbation of immune function may increase susceptibility to infection, impair vaccine responses, and disrupt immune homeostasis. Regulatory guidance, including the ICH S8 Immunotoxicity Guideline, recommends a weight-of-evidence approach in which observations from conventional toxicological endpoints are integrated with functional immune assays to support interpretation of immune system effects. The present study applied an integrated immunotoxicity evaluation framework to examine concordance among structural, functional, and cellular immune endpoints in male Sprague-Dawley rats using a well-characterized immunosuppressive reference compound. Hematological evaluation revealed leukopenia characterized primarily by lymphocyte depletion. Reductions in spleen and thymus weights were accompanied by histopathological evidence of lymphoid depletion in multiple immune tissues, including spleen, thymus, lymph nodes, Peyers patches, and bone marrow. Functional immune competence was assessed through hemagglutination antibody response to sheep red blood cells and delayed-type hypersensitivity assays, both of which demonstrated marked suppression of adaptive immune responses. Flow cytometric immunophenotyping further demonstrated substantial reductions in B-cell populations and decreases in CD4 and CD8 T-cell counts, whereas NK cell populations were comparatively less affected. The concordance of hematological alterations, lymphoid tissue changes, impaired functional immune responses, and lymphocyte subset depletion provides integrated evidence of immune system perturbation. These findings demonstrate that complementary immunotoxicity endpoints collectively support hazard characterization of immune system effects under GLP conditions. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=134 SRC="FIGDIR/small/713556v1_ufig1.gif" ALT="Figure 1"> View larger version (72K): org.highwire.dtl.DTLVardef@beaf9dorg.highwire.dtl.DTLVardef@fb9f10org.highwire.dtl.DTLVardef@187ff06org.highwire.dtl.DTLVardef@1780dc2_HPS_FORMAT_FIGEXP M_FIG C_FIG
Gherardini, J.; Akhundlu, A.; Gompels, M.; Verbinnen, A.; Velasco, S.; Knie, U.; Kassir, R.; Cheret, J.; Paus, R.
Show abstract
We have previously shown that the thyroid hormones triiodothyronine (T3) and thyroxine (T4) prolong anagen, mitigate stem cell apoptosis, and stimulate mitochondrial functions in microdissected human scalp HFs ex vivo. To circumvent the systemic adverse effects of T3/T4, we have asked in the current pilot study whether topically applied T3/T4 retains hair growth-promoting properties. To prove this, we have topically treated healthy full-thickness human scalp skin with T3 (1, 10nM) and T4 (1, 10M) for six days in serum-free organ culture, using an HF-targeting vehicle that contains only FDA-approved ingredients. This showed that, at distinct doses, topical T3 and T4 significantly increased the percentage of HFs in anagen, decreased the percentage of proliferative (Ki-67+) cells in the hair matrix, did not promote melanogenesis (as measured by quantitative Warthin-Starry histochemistry), and significantly increased keratin 15 expression in the bulge. Finally, T3 and T4, at low concentrations, increase the expression of the hair growth promoters IGF-1 and FGF-7. The lower concentration of T3 and both of T4 also significantly increases the number of CD31+ endothelial cells, suggesting a pro-angiogenic effect, which is also important for hair growth promotion. These preliminary results strongly suggest that topically applied thyroid hormones promote hair growth in intact human scalp on multiple levels ex vivo. This invites the intermittent pulse application of topical T3 and T4 as a novel therapeutic intervention for managing hair loss disorders associated with telogen effluvium, such as androgenetic alopecia.
Behr, A.-C.; Vettorazzi, A.; Streel, C.; Mertens, B.; Antonissen, R.; Guerreiro, B.; Ventura, C.; Vilela, R. S.; Novak, M.; Zegura, B.; Reith, F.; Oltmanns, L.; Prisyazhnoy, V.; Suessmuth, R.; Silva, M.; Louro, H.; Marko, D.
Show abstract
Alternaria toxins are naturally occurring food contaminants with limited and often inconsistent genotoxicity and mutagenicity data. Within the European Partnership for the Assessment of Risks from Chemicals (PARC), an OECD-aligned in vitro testing strategy was applied to fill existing data gaps and to characterize the genotoxic potential of major Alternaria toxins using high-purity test materials. Mutagenicity was assessed using bacterial reverse mutation test (OECD TG 471) and SOS/umu assay, while chromosomal damage was assessed using the in vitro micronucleus (MN) assay (OECD TG 487) in TK6 and HepG2 cells, complemented by fluorescence in situ hybridization (FISH) and {gamma}H2AX assay in HepaRG cells. Alternariol (AOH), alternariol monomethyl ether (AME), and altertoxin-I (ATX-I) showed clear mutagenicity in bacteria, whereas altenuene (ALT), tenuazonic acid (TeA), and tentoxin (TEN) were negative under the tested conditions. In mammalian cells, AOH, AME, and ATX-I induced MN formation in TK6 cells at concentrations [≥]5.5 {micro}M, [≥]2.5 {micro}M, and [≥]0.21 {micro}M, respectively, with FISH analysis supporting a clastogenic mode of action. In HepG2 cells, all tested toxins induced chromosomal damage, with effect threshold ranging from [≥]6.25 {micro}M (AOH) to [≥]50 {micro}M (TeA). {gamma}H2AX induction confirmed DNA damage for AOH and ATX-I, and at higher concentrations for TeA (1000 {micro}M). Overall, the data indicate clear in vitro genotoxic potential for AOH, AME, and ATX-I and provide evidence of chromosomal damage for ALT, TEN, and TeA, thereby reducing critical data gaps for hazard assessment.
Borsos, E.; Descamps, B.; Hetzschold, N.; Varga, E.; Marko, D.; Aichinger, G.
Show abstract
The foodborne mycotoxins alternariol (AOH) and alternariol monomethyl ether (AME) have been associated with several adverse effects, including cytotoxicity, genotoxicity, endocrine disruption, and immunomodulation. As these endpoints are typically observed in vitro at micromolar concentrations, the question arises whether such levels are attainable in exposed humans. To address this data gap in chemical risk assessment, a physiologically based kinetic (PBK) model was developed to predict internal exposure doses to AOH and AME in humans. As input parameters, kinetic constants for hepatic glucuronidation were obtained in vitro by incubating Sprague Dawley rat and human liver S9 fractions with 0.5-50 M AOH and 0.5-20 M AME, demonstrating rapid biotransformation in both species. Intestinal absorption of AME and physicochemical parameters were estimated using quantitative structure-activity relationship (QSAR) models. Sensitivity analysis identified parameters describing hepatic glucuronidation and gastrointestinal uptake as among the most influential, confirming the importance of their reliable estimation. The PBK model was evaluated against available rodent toxicokinetic data and subsequently extrapolated to humans. Ultimately, the currently available exposure estimates published by EFSA in 2016 were applied to predict target tissue concentrations, which were compared to points of departure (PoDs) for relevant toxicological endpoints. Even in the most susceptible group of male toddlers, predicted internal concentrations (10-4 M range) were approximately four orders of magnitude below the respective PoDs. Consequently, under the applied exposure assumptions and considering the compounds as isolated chemicals, AOH and AME are not expected to reach systemic or tissue concentrations associated with the investigated effects.
KIM, C. H.; Zhu, Z.; Tagmount, A.; Barbazuk, W. B.; Bacher, R.; Stuchal, L. D.; Martyniuk, C. J.; Vulpe, C. D.
Show abstract
Inter-individual and population variability in susceptibility to chemical exposures confounds determination of threshold exposure levels to protect the most vulnerable. Current risk assessment frameworks, in the absence of empiric chemical-specific data, generally recommend default or probabilistic adjustment factors to account for such variability. We present an experimental approach to incorporate common genetic variants potentially impacting population-level differences in toxicant susceptibility into human cell-based models for any cellular apical endpoint of interest. We focus on the genes with the most common aggregate loss-of-function (LoF) alleles in the gnomAD v3.0 data which we designated as the PopVarLoF set. Unexpectedly, enrichment analysis of these genes found significant overrepresentation of gene products playing important functional roles in toxicology. Interrogation of GWAS and PheWAS databases found that these genes are associated with diverse metabolic phenotypes consistent with the relevance of the PopVarLoF set in studying variability of toxicant response in human populations. We further characterized the PopVarLoF set by developing custom lentiviral CRISPR knockout libraries targeting the PopVarLoF genes to assess their functional essentiality in the HepG2/C3A cell line. Functional disruption of 14 of the PopVarLoF genes ([~]1 %) without toxicant exposure resulted in significant growth defects in this cell line, consistent with the majority of PopVarLoF gene products having non-essential roles. The development of human cell-based toxicity assays or other NAMs which include the empiric assessment of common genetic sources of population variability in susceptibility to chemical exposure could contribute to more robust risk assessment which protects vulnerable populations while reducing uncertainty. Impact statementWe characterize common loss of function genetic variants which could impact toxicant susceptibility and describe an approach to incorporate them into NAMs to enable empiric estimates of the contribution of genetic variability to diverse toxicity endpoints.
Ireland, D.; Zhang, S.; Bochenek, V.; Hsieh, J.-H.; Rabeler, C.; Meyer, Z.; Collins, E.-M.
Show abstract
Organophosphorus pesticides (OPs) are a chemically diverse class of commonly used insecticides. Epidemiological studies suggest that low dose chronic prenatal and infant exposures can lead to life-long neurological damage and behavioral disorders. While inhibition of acetylcholinesterase (AChE) is the shared mechanism of acute OP neurotoxicity, OP-induced developmental neurotoxicity (DNT) can occur independently and/or in the absence of significant AChE inhibition, implying that OPs affect alternative targets. Moreover, different OPs can cause different adverse outcomes, suggesting that different OPs act through different mechanisms. These findings emphasize the importance of comparative studies of OP toxicity. Freshwater planarians are an invertebrate system that uniquely allows for automated, rapid and inexpensive testing of adult and developing organisms in parallel to differentiate neurotoxicity from DNT. Effects found only in regenerating planarians would be indicative of DNT, whereas shared effects may represent neurotoxicity. We leverage this unique feature of planarians to investigate potential differential effects of OPs on the adult and developing brain by performing a comparative screen to test 7 OPs (acephate, chlorpyrifos, dichlorvos, diazinon, malathion, parathion and profenofos) across 10 concentrations in quarter-log steps. Neurotoxicity was evaluated using a wide range of quantitative morphological and behavioral readouts. AChE activity was measured using an Ellman assay. The toxicological profiles of the 7 OPs differed across the OPs and between adult and regenerating planarians. Toxicological profiles were not correlated with levels of AChE inhibition. Twenty-two "mechanistic control compounds" known to target pathways suggested in the literature to be affected by OPs (cholinergic neurotransmission, serotonin neurotransmission, endocannabinoid system, cytoskeleton, adenyl cyclase and oxidative stress) and 2 negative controls were also screened. When compared with the mechanistic control compounds, the phenotypic profiles of the different OPs separated into distinct clusters. The phenotypic profiles of adult vs regenerating planarians exposed to the OPs clustered differently, suggesting some developmental-specific mechanisms. These results further support findings in other systems that OPs cause different adverse outcomes in the (developing) brain and build the foundation for future comparative studies focused on delineating the mechanisms of OP neurotoxicity in planarians.
Quartermain, E.; Zhang, J.; Marczylo, T.; Gant, T. W.; Jacobs, M. N.
Show abstract
Cytochrome P450 (CYP)-mediated biotransformation of endogenous and xenobiotic substances can lead to altered exposure, toxicological impact, or adverse drug reactions. CYP induction data are fundamental to regulatory chemical toxicity hazard assessment because they directly affect the in vivo fate of xenobiotics, potentially influencing their safety and efficacy of pharmaceuticals, and impacting the safety assessment of industrial chemicals, and environmental contaminants. Here we report on the third laboratory supplementary validation of an established and previously validated human HepaRGTM in vitro method able to detect CYP1A2, CYP2B6, and CYP3A4 induction, to support the expansion of the chemical applicability domain beyond pharmaceuticals. This study was conducted to support the part 1 study with additional robust data. We established the test method in-house using the 10 previously validated pharmaceutical proficiency chemicals, then tested a further 6 proposed augmentation chemicals, tebuconazole, benfuracarb, chlorpyrifos, N, N-Diethyl-meta-toluamide, fipronil, permethrin, as tested in part 1, and then four additional chemicals: prochloraz, atrazine, pyrimethanil, and chlorpyrifos-methyl. LC-MS/MS was utilised to measure the conversion of a cocktail mixture of prototypical selective CYP probe substrates to their metabolites, in parallel with mRNA measurements. We achieved high concordance with expected classifications for proficiency and additional chemicals. Comparisons with mRNA-based measurements suggested gene expression may serve as a cost-effective pre-screening tool for CYP1A2 and CYP3A4, though with greater uncertainty for CYP2B6. The data support the robustness of the HepaRG method for CYP induction testing and the adoption of the test method in 2026 as an Organisation for Economic Cooperation and Development Test Guideline. Plain language summaryCytochrome P450 (CYP) enzymes metabolize drugs, pesticides, and other chemicals. Chemicals that increase or decrease CYP enzyme activity can change internal exposure levels, potentially leading to unexpected toxicity or impact drug effectiveness. Reliable in vitro methods to assess CYP induction are needed for regulatory chemical safety assessment. This study describes results from a third laboratory applying a previously validated human HepaRG cell-based method to assess induction of CYP1A2, CYP2B6, and CYP3A4. After successful in-house implementation using ten reference pharmaceutical compounds, the method was extended to ten more industrial chemicals. CYP induction was evaluated by measuring enzyme activity and changes in gene expression. The test method showed a high level of agreement with expected induction outcomes. Gene expression data supported enzyme activity results, particularly for CYP1A2 and CYP3A4. These results strengthen confidence in the robustness and wider applicability of the method for Organisation for Economic Cooperation and Development Test Guideline adoption.
Weller, S. R.; Burnell, J. E.; Aho, B. M.; Obeng, B.; Ledue, E. L.; Shim, J. K.; Hess, S. T.; Gosse, J. A.
Show abstract
People are exposed to high concentrations of antibacterial agent cetylpyridinium chloride (CPC) via personal care and food products, despite little information regarding CPC effects on eukaryotes. CPC is used as an antibacterial agent via a detergent mechanism when above [~]600- 900 M. While three previous studies suggested CPC mitochondrial toxicity, this phenomenon is not well-studied. Here, we show that low-micromolar CPC inhibits mitochondrial ATP production in primary human keratinocytes, mouse NIH-3T3 fibroblasts, and rat RBL-2H3 immune mast cells, in galactose media, which causes cells to produce ATP via mitochondria. ATP inhibition via CPC (EC50 1.7LJ{micro}M) is nearly as potent as that caused by canonical mitotoxicant CCCP (EC50 1.2LJ{micro}M). CPC inhibition of oxygen consumption rate (OCR) tracks with that of ATP: OCR is halved due to 1.75 M CPC in RBL-2H3 cells and 1.25 M in primary human keratinocytes. Here we demonstrate that CPC is more potent as a mitotoxicant than as an immune mast cell signaling inhibitor, an effect published previously. Mitochondrial [Ca2+] changes can cause mitochondrial dysfunction. Here we show, using a novel plate reader assay with reporter CEPIA2mt, that CPC causes mitochondrial Ca2+ efflux from mast cells via an ATP-inhibition mechanism. Using super-resolution microscopy (fluorescence photoactivation localization) in live cells, we have discovered that CPC causes mitochondrial nanostructural defects in fibroblasts, including the formation of spherical structures with donut-like cross section, as quantified by novel Fourier transform analysis. This work reveals CPC as a mitotoxicant despite widespread use, highlighting the importance of further research into its toxicological safety.
Saarimaki, L. A.; Morikka, J.; Pavel, A.; Korpilahde, S.; del Giudice, G.; Federico, A.; Fratello, M.; Serra, A.; Greco, D.
Show abstract
Mechanistic toxicology has emerged as a powerful framework to inform on the safety of chemicals and guide the development of new safe-by-design compounds. Although toxicogenomics provides support towards mechanistic evaluation of chemical exposures, the implementation of toxicogenomics-based evidence in the regulatory setting is still hindered by uncertainties related to the analysis and interpretation of such data. Adverse Outcome Pathways (AOPs) are multi-scale models that link chemical exposures to adverse outcomes through causal cascades of key events (KEs). The use of mechanistic evidence through the AOP framework is actively promoted for the development of new approach methods (NAMs) and to reduce animal experimentation. However, in order to unleash the full potential of AOPs and build confidence into toxicogenomics, robust and unified associations between KEs and patterns of molecular alteration need to be established. Here, we hypothesised that systematic curation of molecular events associated with KEs would enable the modelling of AOPs through gene-level data, creating the much-needed link between toxicogenomics and the systemic mechanisms depicted by the AOPs. This, in turn, introduces novel ways of benefitting from the AOP concept, including predictive models, read-across, and targeted assays, while also reducing the need for multiple testing strategies. Hence, we developed a multi-step strategy to annotate the AOPs relevant to human health risk assessment. We show that our framework successfully highlights relevant adverse outcomes for chemical exposures with strong in vitro and in vivo convergence, supporting chemical grouping and other data-driven approaches. Finally, we defined and experimentally validated a panel of robust AOP-derived in vitro biomarkers for pulmonary fibrosis.
Durnik, R.; Juchelkova, T.; Hecht, H.; Winkelman, L. M. T.; Beltman, J. B.; Comoul, X.; Jornod, F.; Audouze, K.; Blaha, L.; Bajard, L.
Show abstract
As toxicology shifts towards non-animal testing, quantitative models are essential to predict adverse health effects from molecular or cellular perturbations. Quantitative Adverse Outcome Pathways (qAOPs) represent such models, building on mechanistic knowledge and quantifying the Key Event Relationships (KERs) described in AOPs. Despite the recognized need, the number of qAOPs remains limited. Bayesian-based approaches are often chosen for developing qAOP for their flexibility, but most use discretized variables, limiting their predictive power. In addition, these models are mainly built from newly generated data, underexploiting the large amount of information available. This study successfully leverages data from public literature and presents an innovative framework based on continuous variables to develop a Bayesian-based quantitative model for a central KER towards liver fibrosis. The model predicts the probability of the expression fold change for two key markers of hepatic stellate cell activation (aSMA and COL1A1), given the effects on tissue injury, using in vitro data from 9 chemicals. We propose a newly developed workflow to assist in knowledge identification, organization, and extraction from scientific literature and chemical databases. Based on in vitro data and in vivo information from the Open TG-GATEs (Toxicogenomics Project-Genomics Assisted Toxicity Evaluation System) database, we estimate a biologically relevant range in COL1A1 fold change that indicates an activated state of stellate cells and high liver fibrosis odds ratios. Our study provides a case example of integrating published data and continuous variables to build a Bayesian-based model, which constitutes an essential step for predicting liver fibrosis from in vitro data.
Canzler, S.; Lehmann, J.; Schor, J.; Busch, W.; Hackermuller, J.
Show abstract
The grouping of chemicals based on shared properties or molecular mechanisms of action is pivotal for advancing regulatory toxicology, reducing data gaps, and enabling cumulative risk assessments. This study introduces a novel framework usingChemical-Gene-Phenotype-Disease (CGPD) tetramers derived from the Comparative Toxicogenomics Database (CTDbase). Our approach integrates toxicogenomics data to identify and cluster chemicals with similar molecular and phenotypic effects across diverse categories, including pesticides, pharmaceuticals, and industrial chemicals such as bisphenols and per- and poly-fluoroalkyl substances (PFAS). We validated our method by comparing CGPD Tetramer-based clusters with cumulative assessment groups (CAGs) for pesticides, demonstrating strong overlap with established groupings while identifying additional compounds relevant for risk assessment. Key examples include clusters associated with endocrine disruption and metabolic disorders.By bridging omics-derived molecular data with phenotypic and disease endpoints, this framework provides a comprehensive tool for chemical grouping and supports evidence-based regulatory decision-making, facilitating the transition to next-generation risk assessment methodologies.
Chennoufi, M. M.; Dridi, D.; Lasram, K.; Ben Abdeljalil, N.; Omezzine, A.; Boughattas, N. A.
Show abstract
PurposeIfosfamide (IFO) is an effective chimiotherapeutic agent for sarcomas and germ-cell tumors but its clinical use is limited by severe toxicities, particularly hemorrhagic cystitis and encephalopathy. We recently demonstrated that IFO displays a circadian rhythm of lethal toxicity at LD50, with survival, body weight loss and core temperature strongly dependent on dosing time (Chennoufi MM & Boughattas NA, 2025). That study established the presence of an intrinsic chronotolerance rhythm, but lethality prevented organ-specific analyses. Here, we extend this work by assessing whether sublethal IFO exposure (LD30). Methods160 male Swiss albinos mice were synchronized to a 12:12 h light-dark cycle. IFO LD30 was administered at four circadian times: 1, 7, 13, and 19 Hours After Light Onset (HALO). Endpoints included hematology, hepatic enzymes, histopathology (brain, liver, kidney, and bladder) and neurobehavioral function. Organ-specific chronotoxicity patterns were compared to the LD50 chronotolerance rhythm previously reported. ResultsMarked circadian rhythms were observed across all endpoints. Dosing at 7 HALO induced the most severe overall toxicity, while administration at 13 HALO significantly reduced hepatic, renal, and bladder injury, encephalopathy-related lesions appeared only at 19 HALO. Hematological suppression and hepatic enzyme elevations also varied by dosing time, with lowest toxicity consistently at 13 HALO. Neurobehavioral impairment followed the same pattern. The optimal tolerance window overlapped with the LD50 peak tolerance found in our earlier study. ConclusionSublethal IFO toxicity is strongly time-of-day dependent. Administration at 13 HALO (early active phase) minimizes brain and multi-organ injury and aligns with the chronotolerance peak previously defined at LD50. These complementary findings support the development of circadian-based IFO chronotherapy to improve clinical safety.